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ABSTRACT 

A continuous system such as a galactic disc is shown to be well approximated by 
an N-ring differentially rotating self-gravitating system. Lowest order (m=l) non- 
axisymmetric features such as lopsidedness and warps are global in nature and quite 
common in the discs of spiral galaxies. Apparently these two features of the galactic 
discs have been treated like two completely disjoint phenomena. The present analysis 
based on an eigenvalue approach brings out clearly that these two features are fun- 
damentally similar in nature and they are shown to be very Low frequency Normal 
Mode (LNM) oscillations manifested in different symmetry planes of the galactic disc. 
Our analysis also show that these features are actually long-lived oscillating pattern 
of the N-ring self-gravitating system. 
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1 INTRODUCTION 

A typical N-ring differentially rotating self-gravitating sys- 
tem in astrophysical context is the stellar galactic disc or the 
gaseous disc. Disk galaxies often harbour a lot of 'beautiful 
and conspicuous' mostly non-axisymmetric patterns in the 
form of the spectacular spiral structures, bars, rings, warps, 
lopsidedness, scalloping. All these patterns have been a mys- 
tery for a long time and they continue to be even today! It 
is not obvious why these patterns are formed, how they are 
formed and how they are maintained in the discs of spiral 
galaxies? Depending on the orientations of these disc galax- 
ies in the sky with respect to the line-of-sight (LOS) various 
patterns become available for direct astronomical observa- 
tions. In cylindrical polar coordinate system (R, ip, z), suit- 
able for the disc geometry, perturbations describing many 
such patterns are of the type ~ e lmv , ip being the azimuthal 
coordinate and m is the azimuthal wavenumber. Lowest or- 
der modes especially the m=l are of great interest in astro- 
physics. For perturbations normal to the equatorial plane of 
the disc, m=l represents a very commonly observed config- 
uration of the disc known as warps. Warps are systematic 
bending (Lynden-Bell 1965; Hunter & Toomre 1969) of the 
galactic equatorial plane on either side of the galactic cen- 
tre in the form of a characteristic anti-symmetric integral 
sign shape (when viewed with an inclination angle i ~ 90° 
between the LOS and the minor axis of the galactic disc). 
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Whereas in the equatorial plane of the disc (i ~ 0°) such 
m=l perturbations represent another very commonly ob- 
served phenomenon called the lopsidedness. Lopsidedness is 
referred to as the large scale asymmetry in the observed dis- 
tribution of the stellar light (Block et al. 1994; Rix & Zarit- 
sky 1995) as well as in the distribution of atomic hydrogen 
gas (Baldwin et al. 1980; Richter & Sancisi 1994; Hynes et 
al. 1998). Statistical analysis of various observations coming 
from a broad range of the electromagnetic spectrum (visible 
- IR - radio) have established the fact that warps and lop- 
sidedness are ubiquitous in disc galaxies. They are almost 
proving to be an universal pattern in spiral galaxies. Being 
so abundant in nature, such astounding patterns put for- 
ward a question if they are primordial or repeatedly created 
along their evolutionary path. Since the strong differential 
shear in the galactic disc would wash away any such coher- 
ent pattern leading to a phase-mixed equilibrium within a 
billion year time scale, the survival of these patterns intact 
over 10 billion years is a more serious issue. Despite being so 
common there is still no fully convincing and unquestionable 
mechanism for the origin of these patterns in galactic discs. 
Warps and lopsidedness being an outer disc phenomena can 
be used as a potential probe for investigating the nature of 
gravity in the outskirts of the galactic disc, the structure of 
the dark matter halos. These patterns can also put valuable 
constraints on the galaxy evolution and galaxy formation 
processes (Ostriker & Binney 1989; Binney 1992; Combes 
2004). So it is important to have a deeper understanding 
about the nature of these patterns in the galactic disc. Even 
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though warps and lopsidedness are two completely diverse 
phenomena, normally observed in two different class of ob- 
jects face-on (i ~ 0°) and edge-on (i ~ 90°) galaxies, there 
exists an intrinsic similarity between the two observed pat- 
terns. The aim of the present letter is to show that these two 
different phenomena are fundamentally similar in nature. 
We attempt to present an analysis based on an eigenvalue 
approach to bring out the fact that these phenomena are 
basically discrete low frequency normal mode oscillations of 
the galactic discs. 

We consider a razor thin, initially axisymmetric cold self- 
gravitating disc having an exponential distribution of mass 
surface density, E°(7?) = T,oe~ R ^ Rd , with a scale length Rd 
and central surface density Eo. The unit of frequency is 
\J 7rGSo /Rd ~ 53.0 kms _1 kpc _1 for a typical disc like that 
of Milky Way with E = 640.9 M Q pc~ 2 and R d =3.2 kpc. 
Note that these parameters are chosen to represent a typical 
disc model with disc mass, M d ~ 4xl0 10 M©. The disc is 
differentially rotating with angular speed Sl(R) so that the 
unperturbed velocity field is given by v = (0, RSI, 0). In the 
linear approximation for a razor thin disc (z/R << 1) the 
vertical dynamics is almost completely decoupled from the 
planar one and as a consequence we are justified in treat- 
ing lopsidedness and warp as two independent phenomena 
in the disc. 



in order to produce a self-consistent solution of the prob- 
lem. This we achieve using the integral form of the Poisson 
equation for the perturbed disc under the imposed m = 1 
lopsided mode: 
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Where the kernel in eq.(4) is given by 
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The kernel represents the softened self-gravity of the pertur- 
bation, b being the softening parameter. This allows us to 
perform the numerical integration over the disc by removing 
the singularity at R = R' . The second indirect term in the 
kernel arises due to the m = 1 lopsided mass distribution 
about the geometrical centre of the disc. The indirect term 
plays a crucial role in making the disc susceptible to the 
lopsided instability. 

Using the perturbed potential (eq.[4]) and performing 
straightforward mathematical manipulations, we arrive at 
the following integral equation describing the global dynami- 
cal behaviour of the m=l lopsided mode under the collective 
effect of the softened self-gravity: 



2 LOPSIDEDNESS IN A COLD 
SELF-GRAVITATING DISC 

We formulate here the dynamical equations governing the 
global lopsided mode in the razor thin galactic disc by as- 
suming hydrostatic equilibrium along the vertical direction. 
We write the velocity field, surface density and potential of 
the perturbed disc as: vr = v' R , v v = RS~l + v' v , E = E° + E', 
$ = + Then the linearized Euler and continuity equa- 
tions can be written as: 
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In the above equations D/Dt = d/dt + Sid /dtp and k is the 
radial epicyclic frequency in the disc. 

We consider all the perturbed variables as X'(R, tp, t) = 
X'(R)e i(v ^ t) about the equilibrium state. Stability of the 
mode depends on the sign of the imaginary part of to. Lop- 
sided modes in the disc become unstable when Im(oj) > 0. 
They are stable decaying modes when Im(cj)< 0. Modes 
with Im(u;) = are in general stationary van Kampen 
modes; outside the continuum, these modes are pure nor- 
mal mode oscillation of the whole disc. 

Substituting the above form of the perturbed variables 
into eqs.(l-3) and combining the resulting velocity field and 
the perturbed continuity equation in the slow mode limit 
(lu << SI), we get an implicit relation between the perturb- 
ing potential and the perturbed surface density of the disc. 
Now we require the perturbed surface density(E') and the 
potential (<&') to be connected through the Poisson equation 



w 2 E' + S iop (E ')lo + <S i0 p(E') = 



(6) 



Eq.(6) can be solved by recasting it into a matrix-eigenvalue 
problem. By discretizing on a uniform grid with iV radial 
points in the disc resembling a concentric N-ring system we 
can write eq.(6) in a compact form (for details see Saha et 
al. 2007): 
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Where E; op is the eigenvector corresponding to the 
eigenvalue lu. I, Di op , and Si op are the three N x N real 
square matrices. We call I the unit mass matrix, ®i op as the 
general damping matrix and Si op as the stiffness matrix for 
the galactic disc decomposed into an N-ring system. 



3 WARPING OF THE COLD 
SELF-GRAVITATING DISC 

The dynamical equation of a small bending of the disc per- 
pendicular to its unperturbed plane(z = 0) is given by 
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where D/Dt is defined in the previous section and the small 
bending is described by a single function Z\{R, ip, t). Tnaio 
is the vertical restoring force near the disc equatorial plane 
(z — 0) due to a screened isothermal dark matter halo (de 
Zeeuw & Pfenniger 1988). 



Fhaio(R,<P,i) = -v h (R)zi(R,ip,t) 



(9) 
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Figure 1. Complete eigen spectrum of a cold lopsided disc. 
Dashed line denotes lower boundary of the continuum (Q — k). 
The isolated point corresponding to Im(u>) = on the right side 
of the dashed line lies in the P-gap and represent global discrete 
lopsided mode. 
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Figure 3. Complete eigen spectrum of a cold warped disc. 
Dashed line denotes lower boundary of the continuum (SI — v). 
The isolated point corresponding to Im(ui) = on the right side 
of the dashed line is in the P-gap and represents the discrete 
global warping mode of the disc. 
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Figure 2. The pattern speed of a cold lopsided disc. For a disc 
like that of our galaxy, this pattern speed is ~ 0.27 km s _1 kpc -1 . 



where Uh is the local vertical frequency due to the dark mat- 
ter halo. And T se x / is the vertical self- force due to the bent 
disc itself: 
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In the above equations 6 is again the softening parameter. 
zi(R,ip,t) = K{/i u) (i?)e I< ' v_ "*' ) } describes the bending mode 
of the disc. Substituting zi into eq.[8] and discretizing it on 
a uniform grid with N radial points in the disc we have: 
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Where I, D w , and S w are the three NxN real square ma- 
trices. Once again, I is the unit mass matrix, D w is the 
damping matrix and S w is the stiffness matrix (for details 
see Saha & Jog 2006). h w is the eigenvector corresponding 
to the eigenvalue u>. 

Both eq. [7] and eq.[ll] represent a class of nonlinear eigen- 
value problem (here quadratic eigenvalue problem, hence 
QEP) to describe the global behaviour of the m=l modes 
in a self-gravitating disc. 



4 COMPLETE EIGEN-SPECTRUM OF THE 
N-RING SYSTEM 

We treat the galactic disc as a discrete N-ring system 
for studying its two very important global patterns. Both 
eqs.[7,ll] are finite N-dimensional QEP representing a class 
of nonlinear eigenvalue problem. We solve such QEPs by lin- 
earizing it into a 2Nx2N eigensystem of the type Ay = \y. 
In a typical scenario, the disc particles oscillate vertically 
with frequency v and radially with epicyclic frequency k 
while they circulates around the galactic centre with circular 
frequency SI. A local WKB analysis of the m=l type pertur- 
bations in the continuous disc tells us that in the absence 
of any self-gravity any vertical m=l disturbance (warping 
mode) would precess in the disc with two frequencies SI ± v. 
Similarly, any m=l disturbance in the equatorial plane (lop- 
sided mode) of the disc without any self-gravity would pre- 
cess with SI±k. The '+' sign denotes a fast mode while '-' is 
for the slow mode. In galactic dynamics, the slow modes are 
of particular interest because they can survive longer com- 
pare to the fast mode against the differential shear in the 
disc. In normal galactic disc since u, n > SI, the slow mode 
precession frequencies are always negative i.e. they are the 
retrograde modes in the disc. Naturally, according to WKB 
analysis there exists a principal gap (hence P-gap) of length 
2v in the allowed precession frequencies for the warps and 
2k in the case of lopsidedness in the discs of spiral galaxies. 
In absence of self-gravity, there exists no mode in the P-gap 
which is acting like a forbidden gap in the semiconductors! 
The scenario changes qualitatively when there is self-gravity 
associated with each of these patterns in the global sense. 

On solving the N-dimensional QEP (eq.[7]) for the lop- 
sided perturbations, we obtain 2N eigenvalues and these 
eigenvalues are plotted in Fig. 1. We primarily concentrate 
on those slow modes for which Ro(cj) ~ 0. The eigen spec- 
trum shows that there exists a distinct (clearly separated 
from the rest of the points in the Argand diagram) point 
in the P-gap, close to the boundary of the lower continuum 
(SI — k), corresponding to Re(oi) ~ and lm(o>)=0. The 
mode corresponding to this point is an almost non-rotating 
stationary global lopsided pattern of the disc. Being isolated 
in the Argand diagram, this mode resembles like an excited 
normal mode of the disc. Stationary normal mode oscilla- 
tions of spherical stellar system has been studied in detail 
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by many authors (Miller & Smith 1995; Vandervoort 2003) 
previously. The existence of a normal mode of oscillation 
outside the continuum was shown by Mathur (1990) for ID 
and 3D gravitating simplistic system. The presence of such a 
distinct normal mode in our eigen solution for a disc-like 2D 
self-gravitating system confirms the previous trend found in 
theoretical studies namely self-gravitating systems are capa- 
ble of supporting normal modes of oscillation. In Fig. 2, wc 
show that the position of this mode in the complex plane is 
not sensitive to the number of rings used to represent the 
disc once N is moderately larger. It shows that the pattern 
speed (=Re(u>)) of the lopsided mode stays almost constant 
as we increase the number of rings. Thus it would not be 
unfair to conclude that for moderately large N, an N-ring 
discrete system is quite capable of representing the basic 
physics of a continuous system like the galactic disc. Note 
that the above two figures are derived for a purely galac- 
tic disc without any dark matter halo. However, including a 
rigid, non-responsive screened isothermal dark matter halo 
(with Mft ~ 6.3 Md within 10 Rd and producing a reasonable 
flat rotation curve), we find that the spectrum of the lop- 
sided mode does not change qualitatively. As the dark mat- 
ter halo mass increases, the eigen frequency of the isolated 
discrete normal mode (lying in the P-gap) shifts towards the 
lower boundary of the continuum and in the limit of dark 
matter too dominating over the disc, the discrete normal 
mode ceases to exist. But in reality, the dark matter halo 
may not be non-responsive and such live halo might enhance 
the self-gravity of the perturbation which in turn would help 
sustaining the discrete normal mode oscillation lying in the 
P-gap (which we would like to study in future). The spec- 
trum of lopsided N-ring system is quite interesting in the 
sense that it supports a family of unstable modes including 
a very few van Kampen type stationary modes in the con- 
tinuum region. We have checked that the behaviour of the 
eigen spectrum in the continuum region does not change 
with the inclusion of a dark matter halo. 

The complete spectrum of the warped disc is shown in 
Fig. 3. In deriving the eigen spectrum for the warp, we have 
included a dark matter halo of mass ~ 6.3 x calcu- 
lated within 10 R^. The total mass of the galaxy comes out 
to be ~ 3.0 x 10 11 Mq for the assumed disk- halo model which 
produces a flat rotation curve. The real part of the eigen 
values represent the pattern speed of the warping mode. 
The N-ring system produces 2N such pattern speeds. All 
the eigen modes are purely oscillatory. The spectrum con- 
tains a P-gap of length exactly the same size as twice the 
local vertical frequency (2u). Fig. 3 clearly shows the lower 
and upper boundary of the continuum. And there exists one 
very distinct mode in the P-gap resembling an integral sign 
warp as seen in the observations. In the present study, we 
particularly emphasize on the behaviour of the N-ring sys- 
tem representing the actual galactic disc which is basically 
a continuous system. However, even the galactic disc could 
be grainy if we look at it beyond the epicyclic scale which is 
typically much less than the scale length of the underlying 
mass distribution. So unless we go down beyond the epicyclic 
scale the galactic disc more or less can be considered as a 
continuous system. We try to represent the behaviour of 
such a continuous system by an N-ring differentially rotat- 
ing self-gravitating system. 

In Fig. 4 we show how the eigenvalue of the discrete 
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Figure 4. Pattern speed of the warped disc. For a typical disc like 
that of Milky Way, the pattern frequency ~ 0.47 km s -1 kpc -1 . 
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Figure 5. Size of the P-gap in the eigenvalue spectrum of the 
warped disc. 



mode in the P-gap changes as we increase the number of 
rings (N) to represent the continuous disc. In the limit of a 
large N (in statistical sense) the pattern speed of the discrete 
mode does not vary appreciably. The position of the eigen- 
value remains more or less same in the Argand diagram. 
This directly confirms the fact that there exists atleast one 
such discrete normal mode of oscillation of the N-ring sys- 
tem and it stands out in glorious fashion amongst the sea of 
continuum modes to represent a warp of the disc. 
In Fig. 5, we calculate the length of the P-gap as seen in 
Fig. 3. The dots are the one from the eigenvalue analysis. 
As N increases beyond about 100, the gap length evaluated 
from the eigenvalue analysis exactly matches the analytical 
value ~ 0.276 calculated for the assumed disc-halo model. 
This shows a beautiful match between the two. 



5 DISCUSSION AND CONCLUSIONS 

From the present analysis, it is clear that global self-gravity 
plays a crucial role in governing the behaviour of warps and 
lopsidedness in galactic discs. Self-gravity also plays a key 
role in another interesting galactic system known as polar 
rings (PRs) found around SO galaxies. PRs are rotating self- 
gravitating system and differential precession should destroy 
the flatness of the PRs (Sparke 1986; Anraboldi & Sparke 
1994) within a billion year time scale. The fact that we see 
such beautiful polar rings having almost the same inclination 
about the equatorial plane of the host galaxy, implies that 
self-gravity must be playing a major role in maintaining such 
feature. This very basic idea is common in the present anal- 
ysis of warps, lopsidedness and that of polar rings by Sparke 
(1986). So implementing this idea, one can in principle con- 
struct an integro-differential equation for the inclinations of 
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the polar rings/disk and seek a normal mode like solution. 
Basically, understanding of these phenomena then reduces 
to answering the following question: Can we find a bound 
state (here, a discrete normal mode) of the problem in the 
presence of global self-gravity (here, acting like a potential 
well)? Remember, in the absence of a self- gravity, these fea- 
tures are just the free plane waves (lying in the continuum) 
which are going to get wind up due to the strong differential 
precession. 

One of the possible ways to solve such an integral 
equation is to convert it into a matrix eigenvalue problem 
(Kalnajs 1977; Polyachenko 2005). For an infinitely large 
continuous physical system this matrix eigenvalue problem is 
also infinite dimensional. In order to numerically solve such 
an infinite dimensional matrix eigenvalue problem one has 
to truncate it artificially and make it a finite N-dimensional 
matrix eigenvalue problem solvable by computer. Such a 
truncation is meaningful only when the kernel generating 
the matrix is convergent as N becomes larger i.e. the norm 
of the kernel remains finite. So the basic idea here is to see if 
such a finite N-dimensional eigenvalue problem can recover 
the basic physics of the parent continuous system. Our nu- 
merical analysis shows that an N-ring differentially rotating 
self-gravitating system is an adequate representation of a 
continuous system like galactic discs. Note that this idea of 
N-ring has also been used previously in the studies of polar 
rings (Sparke 1986). Apparently, there seems to be a dif- 
ference between the two analyses. Our N-ring system is a 
mathematical construction, while the N-rings (each having 
mass mi) used in Sparke's analysis is physical. In any case, 
the methods should converge as N becomes very large. In 
the linearized version of the two QEPs governing the dy- 
namics of lopsided patterns and warps, the norms of the 
2N-dimensional matrix operators are finite and thus are of 
Hilbert-Schmidt type operators having a discrete spectrum 
of size 2N. Numerical analysis of the eigen spectra shows 
that warps and lopsidedness in a cold disc are fundamen- 
tally similar in nature. Both of them are very slowly (pattern 
speeds are ~ 2 orders of magnitude less than the oribital fre- 
quency) precessing (hence long-lived) normal modes lying in 
the principal gap of the N-ring self-gravitating disc. Such a 
unified picture of warps and lopsidedness as a low frequency 
normal mode oscillation of the N-ring differentially rotating 
self-gravitating system helps us unraveling the more compli- 
cated nature of spiral galaxies. 
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